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Abstract 



Quadrupole excitations of neutron-rich nuclei are analyzed by using the 
linear response method in the Quasiparticle Random Phase Approximation 
(QRPA). The QRPA response is derived starting from the time-dependent 
Hartree-Fock-Bogoliubov (HFB) equations. The residual interaction between 
the quasiparticles is determined consistently from the two-body force used in 
the HFB equations, and the continuum coupling is treated exactly. Calcula- 
tions are done for the neutron-rich oxygen isotopes. It is found that pairing 
correlations affect the low-lying states, and that a full treatment of the con- 
tinuum can change the structure of the states in the giant resonance region. 
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I. INTRODUCTION 



The collective excitations of atomic nuclei in the presence of pairing correlations is usually 
described in the quasiparticle- Random Phase Approximation (QRPA) Although the 
QRPA was applied to nuclear physics more than 40 years ago |§-§J, recently there is a 
renewed interest on its grounds, generated mainly by the studies of unstable nuclei close 
to the drip line. In these nuclei characterized by a small nucleon separation energy, the 
excited states are strongly influenced by the coupling with the quasiparticle (qp) continuum 
configurations. Among the configurations of particular interest are the two-qp states in 
which one or both quasiparticles are in the continuum. In order to describe such excited 
states within QRPA one needs a proper treatment of the continuum coupling, which is 
missing in the usual QRPA calculations based on a discrete qp spectrum. 

In nuclei close to the drip lines one expects also a strong connection between the ex- 
citations of the system and the properties of the ground state, which may present such 
specificities as neutron skins. Therefore, in addition to the qp spectrum, the residual inter- 
action used in QRPA should be determined from the same two-body force as it is done in 
the self-consistent continuum RPA calculations 

In the past years several attempts [jSHToH have been made to describe consistently both 
the pairing correlations and the continuum coupling within QRPA. Thus, in Ref. || a QRPA 
approach was recently developed in which the effect of the continuum is calculated exactly 
for the particle-hole excitations whereas in the particle-particle channel the active space is 
limited to the bound states close to the Fermi level. 

A continuum qp linear response approach in which the continuum is included also in the 



particle-particle channel was studied in Ref. [10], but in the calculations the ground state 



mean field is fixed independently of the residual interaction. 

In this paper we present the first continuum QRPA calculations with the single-particle 
spectrum and the residual interaction determined from the same effective two-body force. 
The ground state is calculated using the continuum HFB approach |]TT[ with the mean 



field and the pairing field described by a Skyrme interaction and a density dependent delta 
force, respectively. Based on the same HFB energy functional we derive the QRPA response 
function in coordinate space. The QRPA response is constructed by using real energy 
solutions for the continuum HFB spectrum. The calculations are done for the neutron-rich 
oxygen isotopes. 

In Section II we present the continuum QRPA formalism, we specialize the corresponding 
equations to systems with spherical symmetry and we discuss the energy-weighted sum rule 
in QRPA. The application of the present theory to neutron-rich oxygen isotopes is done in 
Section III. Section IV contains the concluding remarks. 

II. FORMALISM 

A. Derivation of the generalized Bethe-Salpeter equation 

The coordinate space formalism is naturally adapted to treat properly the coupling to 
the continuum states. In this section we derive the QRPA equations in coordinate space as 
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the small amplitude limit of the perturbed time-dependent HFB equations. We start from 
the time-dependent HFB (TDHFB) equations 

ih™ = [H(t) + r(t),n(t)] (2.i) 

where TZ and TC are the time- dependent generalized density and HFB hamiltonian. The 
external periodic field T is given by 

T = Fe~ iujt + h.c. (2.2) 

where F includes both particle-hole and two-particle transfer operators: 

f = t, F lH c i + E(4 2c H + F ?h c i) ( 2 - 3 ) 

ij ij 

and the particle creation and annihilation operators, respectively. Assuming that 

the external field induces small oscillations around the stationary solution of the HFB equa- 
tions, 

K{t) = 1Z° + H'e- iu}t + h.c. (2.4) 



H{t) =H° + H'e- iuJt + h.c. (2.5) 
the TDHFB equation ( f2.1| ) becomes 

hull' = [W, TZ°] + [H°, W] + [F, K°] (2.6) 
The generalized density variation has the form: 

^ = f R> -n> ) ^ 

where p\, = (0|cjcj|'\ is the variation of the particle density, /€•• = (0|cjCj|') and R\, = 



0\c\c\\') are the fluctuations of the pairing tensor associated to the pairing vibrations and 



3 1 l 

|') denotes the change of the ground state wavefunction |0 > due to the external field. 
Instead of the variation of one quantity in RPA (//) , we now have to know the variations 
of three independent quantities in QRPA, namely p', k' and R'. 

It is convenient to solve Eq. ( |2.6|) in the qp representation in which both 7i° and 7Z° 
are diagonal |J. We have now to express all quantities of Eq. ( p.6|) in this representation. 
The matrix TZ' becomes off-diagonal because of the TDHFB condition TZ' 2 = 71' imposed on 
Eq.(U): 

TZ' ~( ° ^W/ ° v (OR) 

where fa, fa are respectively the qp creation and annihilation operators of an HFB state i 
with energy E^. Consequently, Eq. fl2.6|) gives 
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H' 12 + F- 



12 



K' 12 = 
l] huj - (Ei + E 3 



(2.9) 



K 



in 



H' 21 + F 21 



Huj + (Ei + Ej) ' 

Here, 7i' and F stand for 7i' and F' in the qp representation 

I 12 - 

densities in coordinate representation are denned by: 

p' (rcr) = (0|^ (rcr) ip (rcr) |' 



(2.10) 



We now proceed to calculate H'l 2 and fif^- The variations of the particle and pairing 

are defined by: 

(2.11) 
(2.12) 
(2.13) 



k' (rcr) = (0|^(ra)^(nj) |') 
0|^ f (rcr) ^ (rcr) |' 



K rcr 



where ij}* (rcr) is the particle creation operator in coordinate space and ip* ( r ^) = 
—2ai/;> (r — er) is its time reversed counterpart. The relation between the ip^,ip an d P 
operators is: 



^ (rcr) =^U„ (rcr) fa + V h * (rcr) $ 



(2.14) 



where and are the two components of the HFB wave function of the qp state with 

energy E k . 

Introducing Eq. Q2.14j ) and its hermitian conjugate into Eqs. (|2.11| - |2.13|) one gets with the 
help of Eq. Q) : 



p' a (rcr) = M Tlf + U« 2 (rcr) Tlf, a = 1, 2, 3 



(2.15) 



where we have introduced the following notation for the density variations: 



(p'A 


fp' 


A = 




U/ 





and the 3 by 2 matrices Uij are defined by: 



U t (ra)V 3 (ra) U* (va)V* (rcr) 
UiWUjira) V*(va)V*(va) 
-ViWVjira) -U*(ra)U*(ra) 



(2.16) 



(2.17) 



Here, we have used the same notation as introduced before for the time reversed particle 
operators, i.e., /(rcr)=— 2cr/(r — a). 
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Next, we must calculate the variation Ti! of the HFB hamiltonian in the qp representation. 
This is obtained from the corresponding quantity in coordinate representation through the 
transformation 

H' = W ] H'W (2.18) 

where the matrix W is defined by: 

w =(vu-) < 2 ' 19 ' 

One thus gets: 

W3 3 = / dT E U T (r<r)H nl (r) - Uf (ra)H' 12 (r) - U^\va)H^{v) (2.20) 



(2.21) 



where Uij = —Uji. This antisymmetric combination appears by taking into account the 
relation H! 22 =-H!jl Similar equations hold for the matrix elements of the external field 

(PD- 

In coordinate representation the variation of the HFB hamiltonian is expressed in terms 
of the second derivatives of the HFB energy functional S[p, k, k] with respect to the densities. 
Thus, in our matrix notation we can write (see appendix) : 




(2.22) 



where V is the residual interaction matrix, namely : 

d 2 S 



V Q/3 (ra, rV 



dp/3(r'a')dpa(ra) 



a, 0= 1,2,3. 



(2.23) 



Here, the notation a means that whenever a is 2 or 3 then a is 3 or 2. 

It should be noted that in the three dimensional space, the first dimension represents 
the particle-hole (ph) subspace, the second the particle-particle (pp) one, and the third the 
hole-hole (hh) one. This is due to the definitions ( |2.11| - ^T^ ). 

Using Eqs. (P), ( EICD , (|2~T5| ) and (|2~20|^22|) , we finally land on the coupled equations 



pi = G Vp' + G F 

where F is the three dimensional column vector 




^2.24) 



(2.25) 
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and Go is the unperturbed Green function defined by: 

G ^, r y» = £ W»f(rv) _ ( , 26) 

By definition the QRPA Green function G relates the perturbing external field to the 
density change: 

p' = GF . (2.27) 



Combining with Eq.( 2.24 ) we obtain the generalized Bethe-Salpeter equation : 

G = (1 - GoV)- 1 Go = Go + G VG (2.28) 

In the case of transitions from the ground state to excited states within the same nucleus, 
only the (ph,ph) component of G is acting. If the interaction does not depend on spin 
variables the strength function is thus given by : 

S(u) = --Im [ F n *(r)G n (r, r'; u)F u (r')dr dr' (2.29) 

7T J 

In the equations above we have not introduced explicitly the isospin degree of freedom. 
This can be done directly on the final equations by doubling the dimension of the matrices 
in order to distinguish between neutrons and protons. 



B. Spherical symmetry 

In the case of spherical symmetry each qp state is denoted by the quantum numbers 
(E,l,j,m), where E is the qp energy and (Ijm) are the standard notations for the orbital 
and total angular momenta. Performing the summation over the projection of the total 
angular momentum and over the spin variables one gets for the unperturbed Green function 

G ^(r,r»= £ \A^/g^{ry ]U )Y LM (r)nM(n (2-30) 

LM 

where the expression for the geometrical coefficients ;l j is : 

This expression holds for a residual interaction without gradient terms. We do not give 
here the expressions for the case involving gradient terms since in the calculations of Section 
III we use a Landau-Migdal form for the residual interaction. 

The radial Green functions g?fj p -i q j q { r ^ r '] ^) are expressed in terms of the qp energies Ek 
and the corresponding radial HFB wave functions, i.e., «&(r) = Ui h j h (E k ,r) and Vf.(r) = 
Vi ht j (E k , r). Thus, the radial Green function for a given pair of quantum numbers (l q j q , l p j p ) 
is given by: 
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where 



ft^ - (E p + E q ) +irj hu + (E p + E q ) + irj 



u p (r)v q (r) u q (r)v p (r) 
Upq( r ) = I -u p (r)u q (r) v p (r)v q (r) \ , U pq = U pq + U qp (2.33) 
v p (r)v q (r) -u p {r)u q {r) 



In Eq. fl2.32| ) the J] symbol indicates that the summation is taken both over the discrete 
and the continuum qp states. 

The unperturbed Green function expressed by Eq. (|2.30) can be used together with an 



interaction which does not depend on spin variables in order to get the radial QRPA Green 
function, given by Eq. (|2.28|) . 

All the equations should be written in neutron-proton formalism. Thus, each supermatrix 
(Go, G, V) is divided in 9 blocks corresponding to the ph, pp, hh case and each one of 
these blocks is divided in 4 sub-blocks corresponding to the nn, np, pn and pp quantities. 
The way to calculate explicitly the residual interaction supermatrix is given in appendix. 



C. The energy- weighted sum rule in the QRPA 

It is often stated that the Thouless theorem [Oj concerning the energy-weighted sum 
rule (EWSR) of RPA is also valid for the QRPA. We give here an explicit proof of this 
theorem for the non-trivial case of QRPA. 

The Thouless theorem extended to the QRPA means that the equality : 

Y,K I (u I F | QRPA) \ 2 = \{HFB \ [F, [H, F]} | HFB) (2.34) 

must be satisfied. In the above equation, | HFB) is the HFB ground state and | QRPA) is 
the correlated QRPA ground state while | v) stands for the excited QRPA states. 

To demonstrate this equality we use the QRPA equations in configuration space 0. The 
specific point of the demonstration is that the one-body operator F is now expressed in 
terms of qp operators : 

F = Y. hk'cW = £ fa [U&UwPlPv + U* kl V m ,(5}(5}, + + V^V^p}] (2.35) 

kk' kk' 
W 

The calculation of the l.h.s. of Eq. ( |2.34j ) introduces quantities such as (QRPA | F | v). 
They can be written as : 

(QRPA | F | v) = ( F P )(y») (2 ' 36) 



where 



fw — /, fkk'hkknv, (2.37) 

kk' 
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hkk'w — Vki'Uk'i — VkiUk'v 



(2.38) 



Eq. ( p.36|) allows us to get for the l.h.s. of Eq. ( p. 34 ) an expression similar to the one of 
Ref. p| (Eq. (36)), i.e. : 



(2.39) 



The r.h.s. of Eq. (|2.34|) is obtained by using Eq. ( 2.35 ) and the definition of the A and 
B matrices of QRPA : 



A kk , u , = (HFB\ /3 k ,/3 k , H,/3j/3} t ] \ HFB) 



(2.40) 



B km v = -(HFB | \p v p k , [H,p v ffi\\ \ HFB) (2.41) 

with k < k' and I < I'. Expressing the r.h.s. with A, B and / leads to Eq. (|2.39|) . 

This proof is valid for density-independent force. In the case of a density-dependent force 
the relation should remain valid provided the density-dependent terms commute with F, as 
it has been shown for the RPA case 0. 



III. HFB+QRPA CALCULATIONS OF OXYGEN ISOTOPES 



A. HFB calculations 



We apply our formalism to the calculation of neutron-rich oxygen isotopes 18 > 20 < 22 > 24 O. The 
ground states are calculated within the continuum HFB approach [0] where the continuum 
is treated exactly. The HFB equations are solved in coordinate space with a step of 0.25 
fm for the radial coordinate. In the HFB the mean field quantities are calculated by using 
the Skyrme interaction SLy4 while for the pairing interaction we take a zero-range 
density-dependent interaction given by: 



V, 



pair 



V 



' p{r) 
, Po 



<5 (ri - r 2 ) 



(3.1) 



where Vo, po and a are the parameters of the force. Due to its zero-range this force should 
be used in the HFB calculations with a cutoff in qp energy. To minimize the number of free 
parameters, we adapt here the prescription of Refs. fl4] , |T5"|| which relates the energy cutoff 
with the Vo value for the free neutron-neutron system. To extend this prescription to finite 
nuclei, we use the relation between the energy e of the particle inside nucleus and the energy 
Eq of a free particle: 



e(k)=e u (k)— + U HF 
m* 



(3.2) 



where m* is the effective mass, k the momentum and Uhf the Hartree-Fock potential. Since 
m* depends on the density we take m* /m—0.7 which is the bulk value for m*. From e cuto // 
we can deduce the qp cutoff energy E cuto ff. With this prescription we verified that the 
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calculated HFB neutron pairing gap A n remains constant for each couple (Vo, E cuto ff). In 
the HFB calculations we choose a qp cutoff energy equal to 50 MeV. Then, the prescription 
of Ref. |[3| gives V = -415.73 MeV.fm 3 . The parameter p is set to the usual saturation 



density, 0.16 fm -3 . The value of the parameter a is chosen so as to reproduce the trend of the 
experimental gap. Note that the calculated gap is defined as the integral of the pairing field 
whereas the experimental gap is related to mass differences of the neighboring nuclei and 
therefore, there is no need to have an exact quantitative agreement. We find that the best 
choice is a— 1.5. Note that the trend of the experimental gap is at variance with the empirical 
rule A = 12/\/~A MeV. All the pairing gap values are displayed in Table |. It should be 
noted that, in the case of HFB+Skyrme calculations, if one has a reasonable pairing gap in 
i8,20q th en the pairing is weaker in 22 and absent in 24 0. This is due to the lds/2 and 2si/2 
subshell closure. A similar trend is observed in calculations using Gogny interaction ||16 



As seen in Table I, in the case of QRPA calculations using a Woods-Saxon potential [10 



the gaps are larger due to the fact that the energy distance between the relevant subshells 
is smaller. 



B. QRPA calculations 

In the QRPA calculations the residual interaction is derived in principle from the inter- 
action used in the HFB, i.e., the Skyrme force and the pairing force (|3.1|) . The zero-range 
part of the forces pose no problem. The velocity-dependent terms of the Skyrme force bring 
additional complications which can be avoided by approximating the residual interaction 



in the (ph,ph) subspace by its Landau-Migdal limit [17] where the interacting particle and 
hole have the Fermi momentum and the transferred momentum is zero. The Skyrme inter- 
action has only I = and / = 1 Landau parameters. Taking the Landau-Migdal form for 
the (ph,ph) interaction simplifies greatly the numerical task, at the cost of the loss of some 
consistency. In this work we calculate only natural parity (non spin-flip) excitations and we 
drop the spin-spin part of the (ph,ph) interaction which plays a minor role. The Coulomb 
and spin-orbit residual interactions are also dropped. 

The QRPA Green function can be evaluated starting with the unperturbed Green func- 



tion given by Eq. (|2.32| ). The latter is constructed by using the solutions of the HFB 
equations, i.e., the qp energies and the corresponding wave functions U and V. All the 
qp states are included until an energy cutoff of 50 MeV, allowing pairs of qp energy until 
100 MeV. In a schematic picture, these pairs are representative of excitations from the 16 O 
core, and also excitations of the valence neutrons. The generalized Bethe-Salpeter equation 
( j2.28| ) is solved with a step of 0.5 fm and all radial integrals are carried out up to 22.5 
fm. The strength distribution is calculated until ujMax=50 MeV, with a step of 100 keV 
and an averaging width 77= 150 keV. We have studied two variants of calculations, the full 
continuum variant and a box variant where the qp spectrum is discretized by calculating 
the HFB solutions with a box boundary condition, the box radius being 22.5 fm. For 24 
only box calculations have been performed. 

In a fully consistent calculation the spurious center-of-mass state should come out at zero 
energy. Because of the Landau-Migdal form of the interaction adopted here the consistency 
between mean field and residual qp interaction is broken and the spurious state becomes 
imaginary. We cure this defect by renormalizing the residual interaction by a factor a. We 
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find that in all cases the spurious state J 71 " = 1 comes out at zero energy for a=0.80 . We 
have checked that the EWSR are satisfied within 1% to 5%. 



C. Quadrupole excitations in oxygen isotopes 

We calculate quadrupole strength distributions with the operators F = J2i ff^oifi) 
(isoscalar) and F = J2i r i^2o(ri)t z (i) (isovector). All results presented correspond to the 
SLy4 interaction except when stated otherwise. The strength distributions calculated in the 
neutron-rich oxygen isotopes are displayed in Fig. 1. 

One can identify a strong low- lying state and the giant quadrupole resonance (GQR). 
The low-lying state becomes more isospin-admixed as the neutron excess increases. In the 
case of 24 O the strength distribution is similar to that calculated in Ref. |TIJ with a Woods- 



Saxon potential for the mean field although pairing effects are negligible in our calculation 
whereas the gap A of Ref. flTOJ is sizable. The main difference is the position of the first 2 + 



state located at 4.0 MeV here and 5.0 MeV in Ref. [fTOj] . In the other nuclei this low-lying 
state is at lower energies. This is due to the 2si/ 2 subshell closure in 24 0. The HF single- 
particle energies are given in Table 0. The 2si/2 state is more bound in the 24 O nucleus, 
suggesting a stronger subshell closure in this case. The occupation factors of these states 
calculated in HFB are displayed in Table [TT]. The 2si/2 starts to be significantly populated 
in 22 due to the pairing correlations. In the 18 ' 20 O spectra mainly 3 low-lying peaks are 
present. In the unperturbed case they correspond to the (ld5/2,lds/2), (lds/2,2si/2) and 
(ld5/2,ld3/2) two-qp neutron configurations. Their energies are given in Table |T|. In 18 > 20 O 
the configuration (ld 5 /2,ld 3 / 2 ) has a very low strength whereas the (ld 5 /2,ld 5 / 2 ), (ld 5 /2,2si/ 2 ) 
configurations have similar strength. The effect of the residual interaction, in addition to 
admix the configurations, is to lower the energy of the initial (ld5/2,lds/2) peak and to 
increase the strength of the low-lying state (cf. Fig. [I]). 

The effect of the residual interaction in 22 O is displayed in Fig.|^, showing the isoscalar 
strength functions calculated with the unperturbed Green function Go and with the full G. 
The effect of the residual interaction is to gather the strength to generate collective modes 
such as the GQR and the low-lying state. 



The E2 energy and B(E2) value of the first 2 + state are displayed in Table fV[ As noted 
before, the E2 energy in 24 is larger than in other isotopes due to the 2si/2 subshell closure. 
The El energy in 18 is overestimated and that of 22 is underestimated as compared to 
experiment. This shows, as noted previously in QRPA calculations with a constant gap 
( [pj|l ), that the energy prediction of the low-lying modes is a delicate task in RPA-type 
models. The B(E2) values are well reproduced except for the problematic 18 nucleus. 
This discrepancy observed for the B(E2) value of 18 is also found in several shell-model 



calculations fll9| , |20ll and in previous QRPA calculations |18|, showing a limitation of such 
models to study 18 low-lying states. Indeed this high B(E2) anomaly may be due to the 
presence of deformed states in the experimental low-lying spectra of 18 |21j. Moreover the 
observation of many low-lying states in this nucleus can be described by bands characteristics 
of the rotational spectra. It has been suggested that the low-lying states in these nuclei (such 
as 16 ' 17 ' 18 0) can be described as a mixture between highly deformed states and the usual 



shell model states pl] , p2| . This allows to successfully reproduce both the E2 and B(E2) of 



the low-lying states. As stated in ref. |21|| , for heavier oxygen isotopes, the energies of the 
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deformed states become higher, and thus the admixture is smaller. This may explain why 
the calculated B(E2) are in good agreement with the experimental data for 20 > 22 O nuclei. 
The B(E2) in 24 O is predicted smaller than those of lighter isotopes, which supports the 
2si/ 2 subshell closure effect. In order to display the structure of the low-lying sector, the 
calculated energies and B(E2) of the second and third 2 + states are shown in Table [V] for 
the oxygen isotopes. The 2^ and 2% energies are overestimated in the case of 18 O, whereas 
a good agreement is found for the energy of the 2\ state for 20 O. This may support the 
presence of deformed admixtures in the light neutron-rich oxygen isotopes such as l8 0. 

The calculated M n /M p ratios indicate that the neutron are more coherently contributing 
to the excitation when their number is increasing. For example the M n /M p ratio for 24 
is more than twice the N/Z value, indicating a very strong neutron contribution to the 
excitation. The calculated M n /M p ratio is correctly reproducing that of 20 O deduced from 
proton scattering experiments. In the case of 18 the experimental M n /M p is not well 
reproduced. This is linked to the fact that the B(E2) value is not well described by the 
model. The transition densities for the first low-lying 2 + state of the neutron-rich oxygen 
isotopes are displayed on figure |3|. In the case of 22 > 24 the neutron transition density is 
located more on the surface than the proton one, possibly indicating the presence of a 
neutron skin. 

The pairing effects are depicted in Fig. [| where the continuum QRPA calculation is 
compared to a HF+RPA calculation for the 22 O nucleus. The effect of pairing is to shift 
to higher energy the low-lying peak, and to split the second 2 + state into two states with 
smaller strength. There is also some effect in the GQR region. 

In order to investigate the effect of the density dependence of the pairing interaction we 
have also calculated the strength distributions with a density- independent interaction, i.e., 
Po going to infinity in Eq. fl3.1| ). In the HFB calculation, the Vo parameter has been chosen 
to reproduce the experimental gap of 18 0, Vo=-220 MeV.fm 3 (in this case the prescription of 
Ref. ]T3|] is no longer applied). Fig.|5] compares the results in 18 calculated with the density- 
dependent and density-independent interactions. The effect of the density dependence is to 
increase the energies of the 2 + states, and to slightly lower the strength of the low-lying 
states. 

Box discretization calculations have also been performed in order to test the box bound- 
ary condition approximation. The results are shown in Fig.^|for 22 0. One can see that only 
the low-lying state is nearly insensitive whereas the structure of the GQR is more affected 
by the way the continuum is treated. This shows the necessity of the exact continuum 
treatment in order to study the giant resonances in neutron-rich oxygen isotopes. 

Since our full calculations (HFB+QRPA) have only the Skyrme and the pairing interac- 
tion as inputs the results may be used to learn about the different Skyrme parameterizations. 
Figure 5 shows a comparison of results obtained with the SLy4 SGII and SIII |23| 



for the 20 O nucleus. There is no drastic effect depending on the force. The SIII interaction 
shifts some low-lying states to higher energy and increases the strength of the state located 
around 5 MeV. All three interactions produce a splitting of the GQR but the SGII force 
predicts more strength in the lower component of the giant resonance. 
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IV. CONCLUSIONS 



We have derived the Bethe-Salpeter equation for the QRPA from the small amplitude 
limit of the perturbed time dependent HFB equations. This approach ensures the self- 
consistency at the conceptual level between the mean field, the pairing field and the qp 
residual interaction. The QRPA Green function is decomposed into the ph, pp and hh chan- 
nels. The supermatrix representing the residual interaction is determined self-consistently 
from the Skyrme and the pairing interactions used in the HFB calculations. The Thou- 
less theorem concerning the EWSR sum rule is shown to hold in the case of self-consistent 
QRPA. 

As an application we have studied the quadrupole excitations of the neutron-rich oxy- 
gen isotopes using Skyrme-type interactions for the mean field and a zero-range, density- 
dependent interaction for the pairing field. In the numerical study we have approximated 
the ph residual interaction coming from the Skyrme force by its Landau limit. The cou- 
pling to the continuum appears to have a sizable effect on the GQR and a minor effect on 
the low-lying states. This shows the importance of the full continuum treatment in order 
to study giant resonances in neutron-rich nuclei. The low-lying states are sensitive to the 
pairing interaction. The first 2 + state of 18 is not well described, as previously noted with 
other models such as shell-model calculations. Additional investigations on this nucleus are 
called for. The continuum QRPA shows its ability to reproduce the experimental data of 
the first 2 + state for heavier oxygen isotopes, and it predicts a lowering of the B(E2) value 
for the 24 nucleus. A future improvement of the model will be to include fully the velocity- 
dependent terms of the Skyrme interaction in the ph channel. Work along these lines is in 
progress. 

APPENDIX A: HAMILTON! AN PERTURBATION AND RESIDUAL 

INTERACTION 

The HFB energy functional is written as follow: 

£ = ^tmiPml + ^ £ (ig | V | mp)p pq p m i + J Yl ( lm I I P<l) K lm K pq ( A1 ) 
ml mlpq mlpq 

with V the antisymmetrised interaction and V p the antisymmetrised pairing interaction. 
The HFB hamiltonian is : 

w Ha*^.) (a2) 

with 

9£ A d£ 

^ = %? A - = ^ (A3) 

Next we expand the perturbation of the hamiltonian on the densities perturbations and 
get : 
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_ ff£ t 1 d 2 £ , 1 d 2 £ 

13 ~ \ d^[/ M + 2d^ Ku + 2d^ Kkl (A4) 

d 2 £ 1 d 2 £ 

H 'J = ^ d^/ kl + 2d^/ kl (A5) 

d 2 £ 1 d 2 £ 

^ = ^ d^f kl + 2d^ R ' kl (A6) 

nf = -rq (A?) 

To get TC in coordinate space, we use the relationship between the densities : 

^H = E^H^Hfti (as) 



k ircr 



(A9) 



« M = E 0* ( r ^) M 4 ( A1 °) 

ij 

where fa (re) is the nucleon wave function. 

Using Eq. ( |A4| - |A7| ) together with Eq. (| A8| - |A 1 0| ) allows to calculate 7i' in coordinate 
space, and get Eq. (|2.22|) . 
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TABLES 





A Exp (MeV) 


A Phen (MeV) 


A HFB (MeV) 


A ws (MeV) 


18Q 


1.95 


2.83 


1.96 


2.74 


20Q 


1.83 


2.68 


1.85 


3.13 


22q 


1.52 


2.56 


1.04 


3.30 


24 G 


0.49 


2.45 


0.00 


3.39 



TABLE I. 

Neutron pairing gaps : A Exp is the experimental value taken as the odd-even mass difference 
2~4jl , A Phen is using the empirical 12/y/A MeV prescription f25fl , Ahfb is calculated in the 



present work, and A ws are the gap used in ref |10 





18 


20 O 


22 


24 


!d 5 /2 


-6.7 


0.31 


-6.9 


0.62 


-7.2 


0.93 


-7.7 


1. 


2si/ 2 


-4.0 


0.03 


-4.2 


0.08 


-4.6 


0.18 


-4.9 


1. 


ld 3 /2 


0.3 


0.01 


0.3 


0.02 


0.2 


0.01 


0.2 


0. 



TABLE II. 



Ids/2, 2sx/2 and ld3/2 levels in the 18 > 20 > 22 > 24 O nuclei. For each nucleus the left column shows 
the single-particle energies (MeV) calculated with the HF approximation, and the right 
column displays the occupation factors for the single-qp levels calculated with the HFB 
model. 





18 


20Q 


22 


(ld 5/2 ,ld 5/2 ) 


4.52 


4.16 


4.60 


(ld 5/2 ,2s 1/2 ) 


5.72 


4.36 


3.35 


(ld 5/2 ,ld 3/ 2) 


10.39 


9.09 


7.70 



TABLE III. 

Two qp energies (MeV) of the (ld 5 /2,ld 5 /2), (ld 5 /2,2si/ 2 ), and (ld 5 /2,ld 3 / 2 ) configurations 
of the unperturbed strength function for the 18 > 20 > 22 O nuclei. 
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18q 


20 O 


22q 


24 G 


E2 (MeV) 


3.2 /(2.0 a )) 


2.3 /(1.7 a >) 


1.9 /(3.2 fe )) 


4.0 


B(E2) (e 2 fm 4 ) 


14 /(45 ±2 C )) 


22 / (28 ±2 C )) 


22/ (21 ±8 d )) 


9 


{M n /M p ) 2+ 


2.88 /(1.10 ±0.24 e )) 


3.36 / (3.25 ±0.80 e 


3.53 


4.37 



TABLE IV. 

Energy, proton contribution to the reduced transition probabilities B(E2), and ratio of the 
transition matrix elements M n /M p for the first 2 + state in the 18 > 20 > 22 > 24 O nuclei, calculated 
with the present model. Measured E2, B(E2) values and the M n /M p ratios corresponding 
to the experimental data are displayed in brackets, 
a) ref. p§ ; b) ref. [El ; c) ref. p3 ; d) ref. pi ; e) ref. pi 





is 


20Q 


22 (J 


24 G 


2^ : E2 (MeV) 


5.3 /(4.0 a ^) 


4.2 /(4.1 a )) 


6.2 


7.1 


2+ : B(E2) (e 2 /™ 4 ) 


1.0 


0.3 


1.4 


4.0 


2+ : E2 (MeV) 


9.8 /(5.3 a )) 


8.3 /(5.2 a )) 


7.5 


8.1 


2+ : B(E2) (e 2 /m 4 ) 


1.5 


2.5 


2.5 


0.7 



TABLE V. 

Energy and proton contribution to the reduced transition probabilities B(E2) in the 
is, 20,22, 24q nuc i e i ; calculated with the present model for the 2\ (upper lines) and 2% (lower 
lines) states. Measured E2 values corresponding to the experimental data are displayed in 
brackets. 



a) ref. [26 
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FIGURES 



FIG. 1. Isoscalar (solid) and isovector (dashed) quadrupole strength functions calculated in 
continuum-QRPA for the 18 ' 20 > 22 > 24 O isotopes. 

FIG. 2. Isoscalar quadrupole strength function calculated in continuum-QRPA for the 22 
nucleus. The unperturbed strength (dashed line) is also shown. 

FIG. 3. Neutron and proton transition densities of the first 2 + state of 18 > 20 > 22 > 24 o nuclei. 

FIG. 4. Isoscalar strength function calculated in continuum-QRPA (solid line) and HF+RPA 
(dashed line) with box boundary conditions for the 22 nucleus. 

FIG. 5. Isoscalar strength function calculated with a density-independent pairing interaction 
(solid line) and density-dependent pairing interaction (dashed line) with box boundary conditions 
for the 18 O nucleus. 

FIG. 6. Isoscalar strength function calculated in continuum-QRPA (solid line) and with a box 
discretization (dashed line) for the 22 O nucleus. 

FIG. 7. Isoscalar quadrupole strength function calculated in continuum-QRPA for the 20 O 
nucleus with various Skyrme interactions : SLy4 (solid line), SGII (dashed line) and SIII (dotted 
line) . 
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